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ABSTRACT 



The power spectrum is obtained for the Kolmogorov stochasticity parameter map for WMAP's cosmic microwave 
background (CMB) radiation temperature datasets. The interest for CMB Kolmogorov map is that it can carry direct 
information about voids in the matter distribution, so that the correlations in the distribution of voids have to be 
reflected in the power spectrum. Although limited by the angular resolution of the WMAP, this analysis shows the 
possibility of acquiring this crucial information via CMB maps. Even the already obtained behavior, some of which is 
absent in the simulated maps, can influence the development of views on the void correlations at the large-scale web 
formation. 
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1. Introduction 

Correlations in the full sky or large enough surveys con- 
tain clues to the early Universe and its present structure. 
The acoustic peaks of the cosmic microwave background 
(CMB) power spectrum revealed a set of cosmological pa- 
rameters with particular accuracy (|de Bernardis et al 20001 
IKomatsu et al 2008jl . The baryon acoustic oscillations (see 
Percival et al 2009) are crucial for studies of the forma- 
tion of the large-scale structure, including the role of dark 
matter and dark energy. 

Below, we construct the power spectrum of a novel type 
of full sky map, those representing the distribution of the 
Kolmogorov stochasticity parameter of the CMB tempera- 
ture maps. Kolmogorov's parameter is a descriptor for a de- 
gree of randomness ( [Kolmogorov T933| lArnold 2008aP and 
when applied to the CMB temperature datasets results in 
a map (K-map) ( [Gurzadyan et al 2009 1 that has both fea- 
tures resembling the temperature maps, like the outlined 
Galactic disk, but also ones with different contents. The 
Cold Spot (|Cruz et al 2009p . the non-Gaussian structure 
of negative mean temperature, was noticed thanks to the 
excess of the K-parameter with respect to its mean value 
over the sky. Moreover, the behavior of the K-parameter, 
i.e. of the degree of the randomness was increasing towards 
the boundary of the Cold Spot (Gurzadyan and Kocharyan 
2008, 2009). Both features are compatible to the void na- 
ture of the Cold Spot. Other spots and regions have been 
noticed in the K-map, which are studied with other de- 
scriptors as well (Rossmanith et al 2009), and other no- 
ticed non-Gausianities can also be among the applications 
(Gurzadyan et al 2005, 2008). 

If the Kolmogorov CMB map is able to reflect the 
features in the matter distribution, it is therefore natu- 
ral to study the large-scale correlations in such a map. 



along with the above-mentioned small-scale features. We 
used the latest available full sky maps, i.e. those of the 
Wilkinson Microwave Anisotropy Probe (WMAP) of W, Q, 
V-bands, and the foreground cleaning procedure elaborated 
by Tegmark et al. (2003). The power spectra obtained for 
them have common structures that are, however, absent in 
the simulated maps based on the CMB temperature power 
spectrum. This is the first attempt, and more detailed anal- 
ysis of the K-parameter's power spectra can be performed 
when higher resolution CMB maps are available. 



2. Kolmogorov's stochasticity parameter map 

The Kolmogorov map can be constructed by estimat- 
ing the stochasticity parameter for the CMB tempera- 
ture dataset sequence. Kolmogorov's stochasticity parame- 
ter ( [Kolmogorov 1933] lArno ld 20 08a) is defined for the se- 
quence {Ai, X2, ■ ■ ■ , A„} in increasing order. The cumula- 
tive distribution function is F{x) = P{X < x} , and the 
empirical distribution function is defined as 



rO, x<Xi; 
F„{x)^}k/n, Xk<x<Xk+i, fc = 1,2, . . . ,n- 1 ; 
[1, X^<x. 

The stochasticity parameter is 



A„ = V" sup \F„{x) - F{x) 



(2) 



The universality of this definition stems from how for any 
continuous F , the convergence lim„^oo -P{A„ < A} = 



2 
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<i>(A) , where 



+00 



$(A) = J2 e""" ^ ' = 0, A > , (3) 



k=-c 



is uniform, and <i> is independent on F (Kolmogorov 19331. 

More specifically, to obtain the degree of randomness 
for a given sequence, one must compute the Kolmogorov 
stochasticity parameter A,i, and then the estimated 
Kolmogorov's distribution $ will provide information on 
the degree of randomness in the sequence for the A„ within 
the interval of their probable values, i.e. approximately, 
0.4-1.8 (I Arnold 2008bjl . The mean value of A„ given by 
Kolmogorov distribution is 



A,. 



= J A$'(A)dA « 0.875029. 



(4) 



The behavior of A„ and $ for a set of sequences that 
model the CMB as composition of signals, is studied in 
Ghahramanyan et al. (2009). 

The Kolmogorov map obtained based on this definition 
exhibits, as mentioned above, that structures are linked not 
only to those not iced by other descripto rs but also to those 
indicating voids ( |Gurzadyan et al 2009| . 

3. Power spectrum 

Once the Kolmogorov statistic $ is represented on a map, 
then one can define a correlation function on a sphere in 
spherical coordinates, as for the temperature data, 

C{0) =< $(ni)$(n2) >, nana = cos 6*, (5) 

and expand via spherical harmonics, 

$((?,(^)-^a,™r/™(0,</7), (6) 

where the coefficients a/m, cLS usual, are found from 

a/m = J H9,ip)Y;^{9,^) sin edOd^. (7) 

Then 
and 



0(9) = ^J2^2l+l)CiPi{cos( 



Ci —< a*i„^aim > 



(8) 
(9) 



or 



C, = 



21 



1 ' 



(10) 



■m— — l 



However, for our purposes, i.e. when the (f> is averaged 
within certain numbers of pixels with noise, the cross-power 
spectra C;"* of various bands are more efficient than those of 
autocorrelations; i.e., then one may get more cleaner power 
spectrum for correlations, we study the cross power spectra 
for $ 



"21 



+ 1 



i 



(11) 



ni— — l 



where i = 1,...,8 for Ql, Q2, V1,V2, and W1-W4 

bands. 



In our analysis we used the eight of WMAP's maps, 
of W, V, Q-bands, in the usual HEALPix format 
(jGorski et al 2005)l . of the resolution parameter — 512, 
of a total number of pixels Npix = 12?t-s = 3145728. 
For each — 512 temperature map, we constructed 
Kolmogorov's stochasticity parameter map for Hs = 32, 
Npix — 12288, since for the Kolmogorov map one needs 
about 100 temperature pixels. To obtain the Ug = 32 K- 
map from the Ug = 512 CMB map, each (f> pixel is calcu- 
lated from 64 temperature pixels. 

Then, for the HEALPix map of given riside, the maxi- 
mum I in the obtained power spectrum will be 



In 



'Sims 



(12) 



This corresponds to Imax = 96 for Ug = 32 map and l^ax = 
1536 for Us = 512. The procedure for getting $ cross-power 
spectra included: 

1. calculation of z-th aim for each K-map, 

2. obtaining of all possible combinations of cross-power 
spectra, 

3. estimation of the mean and the error bars for the set 
of spectra: 



0.0025 
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Fig. 1. The mean for 15 cross-power spectra for 
Kolmogorov CMB maps for WMAP's 6 frequency bands, 
V1,V2,W1,W2,W3,W4. 



Cu 



,) = < C'l' >, i^j,0<l< 96; 



(13) 



Note that a/„i-s are complex variables, making the corre- 
lation function complex as well. However, since the noise 
differs from map to map, the resulting complex part due 
to noise is vanishing at cross correlations. The calcula- 
tions were performed for aim without a Galactic disk re- 
gion within ±20°, for 6 and 8 K-maps, and we get 15 and 
28 cross-power spectra, respectively, once their mean and 
error bars were obtained. The results are shown in Fig. 1. 
We see that, for the 28 cross-power spectra, the mean is the 
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same as for 15, but the estimated errors are bigger because 
of using the noisier Ql, Q2 maps. 

The mean power spectrum is similar to the CMB 
pseudo-power spectrum discussed in Hinshaw et al. (2003), 
so one may think to use the Peebles weighting method 
([Peebles 1973: Hivon et al 2002p to find the power spec- 
trum with the Galactic disk. However, this causes two 
types of difficulties. First, we do not have enough pix- 
els {ris = 32, Npix = 12288) to calculate the aim up to 
/ = 250, which is needed for calculating the precise weight- 
ing. Second, even if we keep the Galactic disk region where 
we have approximately <i> = 1, it differs very little from 
the situation if we a priori adopt <I> = 1. The reasonable 
solution seems not to use the Galactic region at all and to 
construct the power spectrum only for odd I, which are not 
affected by the Galactic disk cut. 



4. Foreground cleaned <J> map 

We then obtained the power spectrum of <& using the 
foreground cleaning method developed for CMB maps by 
Tegmark et al. (2003) and the linear combination method 
of (Saha et al. 2006, 2008). This is based on the use of a 
linear combination of different maps with weighting of Wj, 
not only depending on i-th map but also on the multipole 
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where w? is 




0.0000 



Fig. 2. The power spectrum for foreground cleaned 
Kolmogorov maps. 



For constructing of a cleaned $ map, we first calculated 
all aim and then calculated the cleaned a;„i using the rela- 
tion 



E 



w, 



1. 



(15) 



Here C; is an 8 x 8 dimensional matrix constructed by all 
possible auto and cross-power spectra from all maps (see 
eq.(lll), so that Cf^ refers to an inverse matrix. 



and 



(clean) 



(14) 



are 8-dimensional unit vector and its transponated vector, 
respectively. For the covariance of this representation see 
(Tegmark et al. 2003). We get all power spectra in Eq 15 
smoothed by Al = 10 interval to avoid the singular Ui 
matrix. For example, we get triplets of different maps from 
different bands Q1,V1,W1 for wl- 

We thus get 16 different triplets. For any triplet, a lin- 
early superposed aimS is constructed. The last step is to 
find all possible cross-power spectra from those linearly su- 
perposed ones, whose initial map components are different. 
For example, {Ql + Vl + Wl)(g){Q2 + V2 + W2) complies to 
this restriction, but {Q1 + V1 + W1)'S) {Q2 + V2 + WI) does 
not, so only 3 maps of eq.( |15[ ) were used. In this way we 
obtain three cross-power spectra from triplets. The mean 
power spectra from these cross-power spectra is shown in 
Fig. 2. 

5. Simulations 

We repeated the estimations described above for simulated 
maps. Four different types of simulations were constructed 
from: 

a. the real maps' ajm-s (T maps), 

b. real maps with added Gaussian noise of the same pa- 
rameters as the noise in WMAP CMB maps (T-fN maps), 

c. Gaussian maps of the distribution parameters T, a 
from WMAP W band real map (G maps), and 

d. Gaussian maps with added Gaussian noise, both from 
the parameters T, a of WMAP W band map (G+N maps). 

For each group of rt^ = 32 $ simulated map we obtain 
the mean cross power spectra as described above. For 20 
different maps one has 190 cross power spectra. Similarly, 
190 cross-power spectra were computed for the Gaussian 
maps generated with the WMAP's a and mean tempera- 
ture and with superimposed noise of WMAP. 

Although the number of the cross spectra for simula- 
tions is more than those we used for calculating the power 
spectra for real K-maps, neither of the resulted spectra 
shows the features found for real maps with 0.6 and 2.7- 
a level for W and foreground cleaned maps, respectively, 
as shown in Fig. 3. Even more important than the cr-level, 
however, seems that the features only appear at cross and 
not at auto correlations, thus indicating that they do not 
come from the noise in the maps. The principal limitation 
in the above analysis is the angular resolution, since <f> re- 
flects the statistical properties of the signal, the efficiency 
of the method will increase with higher resolution data. 

6. Conclusion 

We have obtained the first power spectrum of Kolmogorov 
stochasticity parameter map of CMB temperature data. 
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Fig. 3. The Kolmogorov power spectra for simulated CMB 
temperature (T+N) and Gaussian (G+N) maps with su- 
perimposed WMAP's noise, averaged over 190 cross-power 
spectra each; the smoothed error bars are shown. 

The WMAP W,Q,V-band datasets were used to compute 
the Kolmogorov's CMB maps. The foreground cleaning 
method of Tegmark et al. (2003) was also applied while 
computing the $ maps. The mean for the set of cross- 
correlated maps was computed. They show features, par- 
ticularly at around I = 25, that are absent in the maps 
simulated either for the WMAP's temperature power spec- 
trum parameters or in the Gaussian maps with superim- 
posed noise, i.e. additional effects to those usually included 
in the simulated models. 

Although the accuracy of the present analysis is limited 
by the WMAP's angular resolution and signal-to-noise ra- 
tio, it shows the principal possibility of obtaining such cru- 
cial information from CMB, and even the already obtained 
behaviors can affect the development of scenarios for the 
void correlations at the large-scale structure formation. 

Higher angular resolution maps expected soon at Planck 
and other experiments will enable the finer structure anal- 
ysis of structures in the power spectra of Kolmogorov CMB 
maps. 
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